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Abstract 

Replica field theory for the Ising spin glass in zero magnetic field is 
studied around the upper critical dimension d = 6. A scaling theory of 
the spin glass phase, based on Parisi's ultrametrically organised order 
parameter, is proposed. We argue that this infinite step replica symme- 
try broken (RSB) phase is nonperturbative in the sense that amplitudes 
of scaling forms cannot be expanded in term of the coupling constant 
w 2 . Infrared divergent integrals inevitably appear when we try to com- 
pute amplitudes perturbatively, nevertheless the e-expansion of critical 
exponents seems to be well-behaved. The origin of these problems can 
be traced back to the unusual behaviour of the free propagator having 
two mass scales, the smaller one being proportional to the perturbation 
parameter w 2 and providing a natural infrared cutoff. Keeping the free 
propagator unexpanded makes it possible to avoid producing infrared 
divergent integrals. The role of Ward- identities and the problem of the 
lower critical dimension are also discussed. 

Spin glasses have posed a formidable theoretical challenge ever since 
Edwards and Anderson H proposed their simple looking model containing 
the two basic features of disordered systems: randomness and frustration. 
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Sherrington and Kirkpatrick (SK) Q denned the mean field version of this 
model on a fully-connected lattice where each spin interacts with all the 
others. The highly non-trivial solution of this mean field model by Parisi Q 
is now generally accepted as the correct one. 

There is, however, no consensus on whether the picture of ultrametri- 
cally organised pure states emerging from Parisi's solution survives in finite 
dimensional, short range systems. "Droplet theory" ||L ||, [/J claims the op- 
posite: in d dimensions the equilibrium spin glass phase is unique, apart 
from reflections of the Ising spins, and a transition from the paramagnet to 
the spin glass takes place only in zero magnetic field. A huge amount of nu- 
merical work has been devoted to clarifying this point (see ^ and references 
therein) and, in contrast to droplet theory, a complex phase space structure 
and ultrametricity seem to emerge in four dimensions. The interpretation 
of recent rigorous results of mathematical physics |], 10] is ambiguous, and 
what we can learn from them is that even the definition of some quantities 
like the probability distribution of overlaps, P(q), is a difficult problem, and 
that chaotic size dependence can render the thermodynamic limit meaning- 
less in some cases. We feel that our findings below about the unavoidable 
infrared problems of replica field theory, arising from a blind application of 
perturbation expansion, are somehow related to these phenomena. 

It is clear that analytical methods are important to settle this decade- 
long debate. One way to depart from the SK limit is the 1/d expansion 
of Georges et al [11] supporting the survival of the replica symmetry bro- 
ken (RSB) picture of Parisi in high spatial dimensions. Another important 
question is whether a spin glass transition exists in a nonzero magnetic field. 
Droplet theory gives a negative answer, while the onset of RSB at the de 
Almeida-Thouless (AT) |l2|] line is an important feature of Parisi's mean 
field theory. Bray and Roberts |R| have studied the problem using Wil- 
son's renormalization group in the replica symmetric (RS) high-temperature 
phase, constraining the fluctuating fields into the replicon subspace. They 
did not find any meaningful fixpoint, thus suggesting that there is no AT- 
line around 6 dimensions. We feel, however, that their projected field theory, 
with all the masses other than that of the replicon taken to be infinite, is not 
sufficient to settle the question since even the crossover region including the 
known zero- field fixed-point [14] is impossible to reach in this model with 
only two cubic coupling constants. Taking into account all the eight cubic 
invariants of the permutation group we would have a large enough parameter 
space with eight y? 3 coupling constants and three masses. A renormalisation 
study, similar to that of Bray and Roberts, of this model could resolve this 
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problem (work is in progress in this direction). 

In this paper we study the replica field theory corresponding to the 
spin glass phase of the Edwards- Anderson model in zero magnetic field and 
just below the transition point. We are mostly interested in systems with 
spatial dimensions around the upper critical dimension d u = 6, the case of 
higher dimensions d > 8 and the problem of the lower critical dimension 
di are left to the end. We have made an extensive perturbation analysis 
of the equation of state and found infrared divergent terms at two-loop 
order for d < 6. We will point out that it is the small mass of the bare 
propagator, rather than the zero (Goldstone) modes, that is responsible for 
these singularities: it is proportional to w 2 , the perturbation parameter. 
When physical quantities are expanded in terms of w 2 , the natural infrared 
cutoff of the theory is destroyed. Nevertheless, below six dimensions, the 
full theory has only one mass scale, and we propose a simple scaling theory 
where all the infinities coming from the infrared integrals are absorbed into 
the amplitudes of scaling forms. 

As explained in Ref. fl5| , bare coupling constants higher than cubic order 
must be set to zero when studying the critical region. The Lagrangean can 
be split into a Gaussian and an interaction part: 

£ = £(2) +£ (/) ; 

with 

ct</3,7<<5 p 

where 

(g-Hp)) = p 2 -r , 

(G- l (p)) n = ~wQ aP , (1) 

and all the other components of G _1 are zero, ro and w are the bare mass and 
cubic coupling constant, respectively, while Q a/ 3 is the exact order parameter 
matrix. Greek indices above and in the following stand for replica numbers 
going from 1 to n, and n is set to zero, according to the replica trick, at 
the very end of the calculations. The above Lagrangean is defined as a 
functional of the fields ip^ with zero average values: 

<<>=0 • (2) 
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(Statistical averages (. . .) are calculated with the weight ~ e .) Eq. || ex- 
presses the shift of (p a P , fluctuating now around the exact Q a f3, as usual in 
the description of symmetry broken phases. Q a p enters also C^: 

£ (/) =£ (D +£ (3) 



with 



£ (3) 



N E ff=0 [ r oQ»f3 + w(Q 2 ) af3 

a</3 

1 



W 



E E^M>- a 



7y ^-^ ^-^ ' Pi Y P2 ^-(pi+P2) 

a<P<jp!,p2 



(JV, the number of sites of the d-dimensional hypercubic lattice, is taken to 
infinity in the thermodynamic limit.) 

We concentrate on the determination of Q a p near T c in zero magnetic 



field where the above Lagrangean is taken with ro 



.(«0 



+ r, r < 1. For 



such a cubic field theory, Eq. 
equations: 



can be written as a closed system of two 



rQ a p + w(Q 2 ) a p 



w 



J E Gap,p P (p) + r^ c) Q af3 



-Y aP ,(3) 
(4) 



where we used the shorthand notation J = J 



(27r) d • 



G is the exact prop- 



agator, while the critical value of ro is to be determined from the implicit 
equation: 



,(<0 



lim 

Q^o Q. 



1 



-it; 



a/3 



^ap,/3p 



00 



(5) 



with limQ^o meaning that all the elements of Q go to zero, depends 



on Q and 7q in a complicated manner. By Eq. [|, is a function of the 
coupling constant w, and it can be calculated perturbatively. 

Computing Q a p from Eqs. |3] and || is a difficult problem, as it is for 
any nontrivial field theory. Without assuming an ansatz for the structure 
of Q a /3, it is completely hopeless to find even an approximate solution. We 
can try, however, to follow the traditional way by starting from a mean 
field-like equation whose solution is taken as a zeroth order approximation 



(c) 
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in a systematic perturbative treatment. The most obvious way of denning 
a mean field equation of state from Eq. |3| is to drop the loop term Y a g. Any 
solution found must then be submitted to the test of stability: the mass 
operator T = G~ 1 (p = 0) should not have negative eigenvalues. The illusion 
created by the marginally stable replica symmetric mean field solution is 
dispelled by calculating the first loop correction (this was made in Ref . ) : 
Bray and Moore found that the hitherto zero, so called replicon, eigenvalue 
moves to a negative value, a clear indication of an effective quartic coupling 
(we call it u) generated by the cubic coupling at the first loop level |1| . 

For a meaningful mean field theory, we have to keep the loop term Y a p 
in Eq. || and make a reasonable truncation in Eq. ||. We can avoid instability 
by assuming an infinite step RSB structure for Q a p. It turns out that in 
this case Q a p is highly insensitive to the details of the propagators in the far 
infrared region, while the near infrared propagators which determine Q a p 
are not influenced by the quartic couplings (for a discussion of the behaviour 
of the propagator components in the far and near infrared, corresponding 



to the two mass scales, see Ref. [15]). We have therefore a freedom in 



truncating Eq. ||, we can even completely neglect when computing the 
average, thus taking G = G. Alternatively, we may take into account the 
quartic couplings, u for instance, being generated from the interaction part. 
In any case, the propagators are the same in the near infrared, while only 
the amplitudes are influenced in the far infrared region giving a relatively 
moderate divergence G ap ^ p ~ p~ 3 . The role of the small mass is to provide 
an infrared cutoff, and our mean field theory is certainly meaningful for 
d > 3. 



Using the tables from Ref. 15] for the near infrared propagators, we can 
approximately compute the Parisi order parameter function in the vicinity 
of T c : 

w 

Qa/3 = Q(x) = —x , x = a<T\(3<xi , (6) 
2u 

Q{x) = Qi = —r , x > xx , 
2w 

x\ = -V , (7) 

where Q\ and x\ stand for the plateau value and breakpoint of Q(x), re- 
spectively. The effective quartic coupling is [17, [0| : 



,.4 " 



p 4 (p 2 + r) 2 ^ ^ 
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Q(x) in Eq. || is very much like that of the SK model, where the coupling 
constants w and u are 1. We can learn further details about Q(x) by writing 
a scaling form for Y a p = Y(x). Considering the two mass scales of the mean 
field propagators 

m large ~ r > 

2 ^2 
m small ~ — r = x l r , 



we have in the critical region: 

Y(x) = w / A 4/mf (-,—,-) , (9) 
J p r x\r x\ 

where / m f is a scaling function. The truncation procedure for the propagator 
has an effect on / m f in the momentum range around the small mass (p 2 ~ 
x\r), but the near infrared region (p 2 ~ r) and the power of the far infrared 
limit (~ p _1 ) are not influenced by it. 

The mean field approximation presented above is highly nontrivial, and 
the order parameter function depends in a very complicated manner on the 
perturbation parameter w 2 . The origin of this complicated dependence can 
be understood from the fact that the small mass is proportional to w 2 : 

"Viall a ~ W 10 

giving rise to a u>-dependence of the propagator G. At this point we are 
tempted to compute Q{x) by perturbation expansion, i.e. order by order in 
w 2 . This, however, inevitably leads to infrared divergent terms since the 
infrared cutoff provided by the small mass is then destroyed. For example, 
in the expansion of Q\ we have a third order term like 



wQi ~ w 6 r d ~ 4 J A ^ 



which is infinite for d < 8. Similarly, as we will see later, when using simple- 
minded perturbation expansion for the full theory we may encounter infrared 
divergent infinite terms that can, however, build up an infrared cutoff, thus 
rendering the theory well-behaved above the lower critical dimension. 

Following standard arguments of the theory of critical phenomena, we 
expect a qualitative agreement between full and mean field theory above the 
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upper critical dimension (d > 6). Near T c , Y(x) [x = a Ci/3) in Eq. || can be 
written as: 



y (*) = ^/ A i/(4'^7>-) > 6<d<8 , (11) 

where the scaling function / depends now on the coupling w and, for later 
reference, the mean field exponents rj = 0, v = 1/2 were included. A 
renormalised u is now understood in Eq. fj] for the breakpoint of Eq. [ll], 
preserving its temperature scaling 

xi~ri" 3 , 6<d<8 . (12) 

Going below the upper critical temperature seems to be easy now: the 



exponents rj and v in Eq. 11 are then nontrivial and depend on the dimension 
(they have been computed from the critical theory up to 0(e 3 ), e = 6 — d, in 
Ref. ]14| ). The most important point to clarify is the role of x\ in the scaling 
of the Y(x) term. Eq. [l^ suggests that x\ becomes temperature independent 
in d = 6. We have made a rather difficult calculation of the 0(e) correction 
to the temperature exponent of x\ (see later) and found that it is zero. The 
hypothesis of a constant, nonuniversal, i.e. temperature independent and w, 
A dependent, x\ is strongly suggested by this finding. The simple scaling 
picture restored in d = 6 would then survive even below the upper critical 
dimension and Y(x) would take the form 

Y W = ™j A ^K^,x) , d<6 , (13) 



with the nonuniversal / replacing / in Eq. [11 . 

Provided that our assumption of a temperature independent x\ is cor- 
rect, a simple form for the order parameter function is strongly suggested^]: 

Q(x) = r?Q(x) , (14) 

where Q(x) depends only on w and A and for the critical exponent j3 we 
have the scaling law 

(5 = \y{d-2 + v ) . (15) 

1 This factorization of r 13 from Q(x) was raised by A. P. Young some years ago in private 
communications. 
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We can find full consistency between Eqs. 13 and |14| by assuming the limiting 
form 

1+/3 

f(v,x)~v~— , u>l (16) 

and using the scaling law Eq. |l5[ Changing the integration variable from p 
to p = p/r v : 

Y{x) = r^w -i-^/(p' 2 , x) r 1+fi Y x {x) + r 2/3 y 2 (x) , 

where the first term is the leading contribution coming from the upper limit 
of the integration while the second one is from the bulk. A Q(x) as in Eq. H 
is then obviously consistent with Eq. |3[ 

The above theory for the spin glass phase below d = 6, despite its con- 
sistency and its clear connection to mean field theory, does require some 
analytical support. We have made a detailed perturbative computation of 
the equation of state H and our findings can be summarized in the following 
points: 

1. A Parisi RSB scheme is building up order by order with a polynomial 
Q{x) = ax + cx 3 + . . .. Consequently, f(v, x) of Eq. 13 must be also a 
polynomial in x for any fixed, nonzero v. 

2. We computed several logarithms of r to check the above scaling theory 
and found agreement with calculations of the critical exponents in the 
critical (massless) phase J14| |. The most important result comes from 
the second order (two-loop) coefficient of x 3 log r that provides the 
temperature exponent of the slope a of Q{x) showing that it is equal to 
(3, at least up to 0(e). This is consistent with the scaling form in Eq. [Ll| 
and supports the basic assumption of a temperature-independent x\. 

The exponent (3 is originally defined by 

Qi = Cr? (17) 

and the one-loop xlogr term gives (3 correctly up to 0(e) {(3 = 1 + |). 
Furthermore, using the two-loop x log 2 r results, exponentiation was 
checked and agreement with the known fixed-point value obtained [0] 
(the latter result has been published in Jlf 



3. We know from exact Goldstone-theorem-like considerations |18|] that 
an infinite-step RSB solution of the equation of state, Eqs. |3| and ||, 
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must be marginally stable since the most dangerous replicon eigenval- 
ues are Goldstone-modes. In contrast to the RS case the one-loop 
calculation gives indeed a marginally stable spectrum of the mass op- 
erator. 

4. The above three points support a Parisi-type structure of Q a p, but 
now we have to deal with the problem of the unremovable infrared 
divergences appearing when a blind expansion in term of w 2 is used. 
Up to second order, Y a p can be written as 

Y a p = E af3 + F af3 + 0(w 5 ) , (18) 



where 



E a /3 = w 



f E G ap Mp) + r { o C) Qap , (19) 



In Eqs. ^ and ^ we use G instead of G, where G is defined like G in 
Eq. 0, but with r and Q^l replacing tq and Q a /3, respectively. Q^l is 
the leading term of the order parameter: 

Qaf3 = Q { ^ + 0(w 2 ) . 

The matrix elements of G -1 are now all zeroth order, nevertheless G 
is not a free propagator in the usual sense: since x± ~ w 2 , all the sizes 
of the blocks in the hierarchical construction are also proportional to 
w 2 and thus G can be expanded in w 2 : 

G = G(°) + G« + . . . . 

As it turns out, G*- ** is nothing but the near infrared propagator whose 
components were listed in Ref. (l^ and further terms can be calculated 
by straightforward, although lengthy, perturbation methods (the com- 
ponent G« l ' was needed for the present calculations). Two types of 
infrared divergences with different origin enter the equation of state: 

• In Eq. 20, where I^ 1 ) is the first order contribution to the two- 
point vertex operator T = G _1 whose zero momentum limit is the 
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exact mass operator, the high infrared power of the two G's could 
lead to infinite terms. These are, however, compensated by the 
zero modes of the mass operator T(0), as a direct consequence of 
exact Ward-identities (very much like in the case of the Heisen- 
berg model). These Ward-identities follow from the permutation 
symmetry of the n replicas which becomes a continuous symme- 
try in the infinite step RSB limit [|l^] . Thus F a p is innocent and 
well-behaved even below d = 6. 

• Ward-identities cannot help us to get rid of infrared danger in 
E a p = E{x) of Eq. [TO]. To illustrate the problem, we write down 
our result for x = x\\ 

A 



E{x x ) = wr 2 J { gi gl-2glg 2 ) + 




where 



-\g\ - h gi gl - 28g 2 g 2 + 



+32gi) + 0(w>) , (21) 



9\ = -z-. — and g 2 = -~ 

p z + r p £ 



All but one term in Eq. |2lj give contributions for the computa- 
tion of j3 (log's), checking of exponentiation (log 2 's) and even to 
the amplitude C in Eq. |l7| (simple powers) up to 0(e 2 ). The 
term — 4g§ is, however, infrared divergent, leading to an infinite 
contribution to C at order e 2 . 

The origin of the infrared terms is identical to that in the mean field case: 
2~2 p =£a,p @ap,Pp(p) is P~ 3 in the zero momentum limit, but when expanding 
it in w 2 , high infrared powers are generated order by order. These enter 



the amplitudes of scaling forms like Eqs. [L3], [L6| and 17, suggesting that 
amplitudes cannot be expanded without producing infinite terms. 

The above scaling theory of the equation of state can be extended to the 
propagator G (a first attempt in this direction has been made in |2^, |2l||). 
Indeed, it is plausible (and calculations on the Gaussian level suggest) that 
any component of G has the scaling form p~ ( - 2 ~ v ^ g(p 2 /r 2u , x, . . .) where the 
scaling function g is, of course, different for different components and . . . 
stands for other possible overlaps of the given propagator component. We 
record here the explicit formula for the diagonal element G^ which is the 
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Fourier transform of the overlap of the spin-spin correlation function (for 
the physical meaning of the propagators, see Ref. ||l5||): 



1 p 2 

GTJ(P) = ^5diag(^7,£) , d<6 . (22) 

For x = 0, it is also the most accessible to numerical simulations [||. 

The theory presented above is based on the assumption of an ultra- 
metrically organized Parisi RSB ansatz. Its analytical justification requires 
perturbation methods like the e-expansion to compute critical exponents (rj, 
v and (5) and check exponentiation. Amplitudes, however, are infected by 
infrared divergences, and the resummation of these infinities seems to be a 
hard task. The difficulties in using a perturbative approach for a system 
with a complex phase space structure has been discussed in the context of 
the random manifold problem in Ref. [19]. Numerical methods would be 



valuable to check formulae like Eqs. 14 and 22, and compare the measured 



exponents with those obtained at or above T c . For this purpose, we make 
further assumptions for G^'* in Eq. [2^ when x = 0, namely 

gdiagCwjO) — > constant , v — > oo , 

SdiagM) ~ v~§ , v^O . (23) 



The value of the exponent A has an important role in the determination of 
the lower critical dimension d\ (see below). 

Up to now, we confined our discussion to systems with d < 8. Above 
eight dimensions we must keep the bare quartic coupling u to avoid instabil- 
ity of the RS solution, and all the formulae containing u are still valid with 
u replacing u. Y(x) in Eq. |3| has the mean field form like in Eq. 11 with 



rj = 0, v = 1/2 and x\ = ur/w 2 ~ r. The scaling function / depends now on 
w, u, A and, in high enough dimension, possibly on other higher order bare 
couplings, and on the l.h.s of Eq. ^ terms like ^uQ^p appear. The theory is 
basically nonperturbative in the sense that a perturbation expansion in the 
bare couplings destroys the infrared cutoff at p ~ r. 

The problem of the lower critical dimension di where the above scaling 
theory blows up, can be approached from two different directions: 

• In the case of any critical transition, T c can reach zero temperature in 
low enough dimension which, in a field theoretic representation, man- 

(c) 

ifests itself in a divergent 7q . From the observation that the second 
term Qapr^ in Y a p (see Eq. ||) is to cancel the leading contribution 
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of the first one when p 2 3> r, we can put Eq. [5] into the simple form 
consistent with our scaling picture: 



+o J p 2 ^ \ r 2u I J pt-n+P/v 
where the asymptotic form f(v) ~ v~^^ 2u , v — > oo, has been used. 



This can now be compared with the subleading term of Eq. 16 



re- 



maining after the cancellation mentioned above. Using the scaling law 



of Eq. 15, we arrive at the conclusion that the integral of Eq. [2j is 
convergent provided d > 2 — r]. The lower critical dimension following 
from this argument satisfies the equation 

df ) =2- V {df ) ) . (25) 

In systems with Goldstone modes, the high infrared powers of the 
propagators lead to divergences in local quantities at some d\ , thereby 

signalling the breakdown of the ordered phase. When d^ > d[ , it 

is the lower critical dimension (a situation known from the XY and 

Heisenberg models where di = 2). As it was suggested in Ref. [EQ], we 

o o — 

must consider the most dangerous component whose far infrared 



behaviour can be obtained from Eqs. 22 and 23: 

g?;?(p)~p- (2 -" +a) , p^O . (26) 

From this we have 

4 2) = 2-r ? (4 2) ) + A(d / {2) ) . (27) 
In Ref. |2GJ, after an incomplete perturbative analysis, it was argued 



that A = /9/z/, giving rise to d^p = , leading to a d\ which is 
definitely smaller than three. Numerical estimate of the exponent A 
for d = 3 can be found in Ref. ||, namely A = 2(3 /v — a with a ~ 0.5 
which gives A ~ 0.1. This result should be taken with some care given 
the small size of the system considered. If accepted, then it gives a 
smaller A than the previous one: A = (3 jv 0.3, d = 3, and is thus 
even more consistent with di < 3. (di = 2.5 was suggested by Franz 
et al p2| as the lower critical dimension where replica symmetry is 
restored.) Nevertheless, we cannot exclude the possibility di = 3, and 



more numerical and analytical work [23] would be useful to settle this 
problem. 
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